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Abstract 

The inclination of low-eccentricity orbits is shown to significantly affect the orbital 
parameters, in particular, the Keplerian, nodal precession, and periastron rotation 
frequencies, which are interpreted in terms of observable quantities. For the nodal 
precession and periastron rotation frequencies of low-eccentricity orbits in a Kerr field, 
we derive a Taylor expansion in terms of the Kerr parameter at arbitrary orbital 
inclinations to the black-hole spin axis and at arbitrary radial coordinates. The particle 
radius, energy, and angular momentum in the marginally stable circular orbits are 
calculated as functions of the Kerr parameter j and parameter s in the form of Taylor 
expansions in terms of j to within 0[j 6 ]. By analyzing our numerical results, we give 
compact approximation formulas for the nodal precession frequency of the marginally 
stable circular orbits at various s in the entire range of variation of Kerr parameter. 
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Introduction 

The X-ray quasi-periodic oscillations (QPOs) discovered in low-mass X-ray binaries 
(LMXBs) commonly show a variety of modes. In particular, there are horizontal-branch os- 
cillations at low frequencies i^hbo 1-100 Hz and two peaks at frequencies V\ and v 2 (~1 kHz) 
in the power spectrum (van der Klis 2000). In many Z-type sources, as well as in several atoll 
sources, the frequencies z/hbo exhibit an almost quadratic dependence on v 2 (Stella and Vietri 
1998; Psaltis et al. 1999). For several sources (e.g., for Sco X-l, 4U 1608-52, 4U 1702-43, 
4U 1735-44, XTE J2123-058), the frequency difference v 2 — v\ = Av decreases with in- 
creasing v 2 . At the same time, for many Z-type (GX 5-1, GX 17+2, Cyg X-2, GX 340+0, 
GX 349+2) and atoll (4r 0614+09, 4r 1636-52, 4r 1705-44, Aql X-l) sources, the difference 
Au is constant, which allows it to be associated with the stellar rotation frequency (van der 
Klis 2000). The frequency v 2 is currently identified with the Keplerian rotation frequency of 
clumps of matter near a compact object [y 2 = vk) i n almost all interpretations. 
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Cui et al. (1998) pointed out that the nodal precession of circular orbitsf] slightly inclined 
to the equatorial plane of a black hole could be of importance in interpreting the stable QPOs: 
for several black-hole candidates and microquasars, they found the predicted frequencies of 
nodal precession to agree (if the Kerr parameter is chosen in the range 0.37 — 0.9) with the 
observed QPO frequencies. 

Stella and Vietri (1998) justified the formula for the nodal precession frequency z/ no d of 
orbits slightly inclined to the equatorial plane and identified its even harmonics with the 
frequency ^hbo- The frequency v\ was associated with the periastron rotation frequency 
^per = i^k — of a low-eccentricity orbit (a general-relativity effect). They showed that the 
theoretical dependence v r {y$) for an appropriately chosen mass of the object in the range 
1.8 — 2.2M was in good agreement with the experimental points in the (1/2, Av) plane for 
various sources with evolving frequencies. In the Kerr solution, the frequency v r becomes 
zero at the marginally stable orbit. Therefore, when a radiating clump of matter passes from 
one Keplerian orbit to another^ by increasing its rotation frequency, the frequency difference 
between the peaks Av = v r decreases and even approaches zero at the marginally stable 
orbit. This tendency was observed for some of the sources (see above). 

Morsink and Stella (1999), Stella and Vietri (1998), and Psaltis et al. (1999) showed 
the theoretical dependences z/ nod (z/ K ) to agree with the measured source frequencies in the 
"2 3 ^hbo plane if z/hbo was identified with an appropriate even harmonic of the nodal fre- 
quency. 

The proper rotation of the source (a neutron star or a black hole) remains a free param- 
eter, and it can be chosen by using the stable frequency observed in LMXBs in the X-ray 
band during outbursts. 

Another free parameter is the orbital inclination s to the equatorial plane of a neutron star 
or a black hole, which significantly affects the observed nodal precession, periastron rotation, 
and Keplerian frequencies near the marginally stable orbit. Below, we give an example based 
on our results. If the mass of a black hole or a neutron star is M, then for the Keplerian 
frequency in the marginally stable orbit = 1.2 (2.2M /M) kHz, the frequency of nodal 
precession u nod is, respectively, 123 (2.2M /M) Hz at s = 0; 93.6 (2.2M /M) Hz at s = 10°; 
65.3 (2.2M /M) Hz at s = 30°; 52.6 (2.2M /M) Hz at s = 45°; 41.18 (2.2M /M) Hz at s = 
80°; and 41.11 (2.2M /M) Hz at s = 90°. The nodal precession frequency changes by almost 
a factor of 3 as the inclination changes from 90° to 0! Thus, by simultaneously measuring 
three quantities, M, z/hbo 3 an d ^2, we can determine the inclination of the marginally stable 
orbit to the spin axis! 

In general, the tilt of an accretion disk to the equatorial plane of a compact object can 
be finite. The Bardeen-Petterson (1975) hypothesis of the accretion-disk transition into the 
equatorial plane breaks down when the radiative forces that twist the inner edge of the disk 
are taken into account. Pringle (1996) showed that a flat disk was unstable to disturbances 
in the presence of a central radiation source. 

Below, we derive formulas for the Keplerian, nodal precession, and periastron rotation fre- 
quencies of low-eccentricity orbits at an arbitrary (finite) orbital inclination to the equatorial 
plane in the form of Taylor expansions in terms of the Kerr parameter. In contrast to previous 
studies, in which calculations were performed for selected values of the constant Q [with an 
unclear physical meaning, as was noted by de Felice (1980), and coinciding with the square 
of the total angular momentum only in the weak-field approximation Q), we consistently 

^ense-Thirring (1918) effect. 

2 Through turbulent viscosity and radiative deceleration. 



use the smallest latitudinal angle 9_ = s reached on a bounded trajectory^ as a trajectory 
parameter. At s = n/2, the derived formulas (51) transform into Taylor expansions of the 
formulas by Okazaki et al. (1987). An analysis of our numerical results has yielded compact 
approximation formulas for the nodal precession frequency of the marginally stable circular 
orbits in the entire range of Kerr parameters for several finite orbital inclinations. 

In contrast to rotating black holes, an intrinsic quadrupole component appears in the 
external fields of rapidly rotating neutron stars (NSs). For NSs with a stiff equation of state, 
this component can be several times larger than the Kerr quadrupole moment. The effect 
of a non-Kerr NS field induced by rapid NS rotation on the parameters of the marginally 
stable orbit and energy release in the equatorial boundary layer was analyzed by Sibgatullin 
and Sunyaev (1998, 2000a, 2000b) using an exact solution for a rotating configuration with a 
quadrupole moment (Manko et al. 1994) [approximate approaches with multipole expansions 
of the metric coefficients at large radii were developed by Laarakkers and Poisson (1998) 
and Shibata and Sasaki (1998)]. Our results refer to nonequatorial, nearly circular orbits 
in the fields of black holes and NSs with a soft equation of state at a moderate rotation 
frequency. Markovic (2000) considered finite-eccentricity orbits in a Kerr field and in the 
post-Newtonian approximation for a field with a finite quadrupole moment. The precession 
of orbits slightly inclined to the equatorial plane in the fields of rotating NSs was numerically 
calculated by Morsink and Stella (1999) and Stella et al. (1999). 



Geodesies in the Kerr solution 

It is well known that the equations of geodesies in Riemannian spaces = can be written in 
Hamiltonian form. The corresponding Hamilton-Jacobi equation is 

Q={g^S ii S ij +l)/2 = 0, S n =— , (1) 
... = 0,1,2,3. 

The generalized momenta Pi = S,i are related to the 4- velocity components by the Hamilton 
equations 

dx j dQ jk dpj dQ dg kl 

1~ = TT~ = 5 ? Pk, ~r = = -jr-PkPi- (2) 

dr opj dr dx d ox 3 

For the Kerr solution in Boyer-Lindquist coordinates, the contravariant metric components 



arc 
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p 2 c 2 A 
1 - 2rGM/c 2 p 2 ) 
Asin6» 



o</> = 2ar tt = ee J_ 

9 ~ p 2 A ' 9 p 2 ' 9 " p 2 ' 

A = r 2 + a 2 - 2rGM/c 2 ; p 2 = r 2 + a 2 cos 2 6. 



3 Bclow, s is called the inclination angle between the angular velocity vector of a compact object and the 
orbital surface for short. 
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The constant a is related to the angular momentum of a rotating black hole described by 
the Kerr solution by the equality J = Mac = GM 2 j/c, where j is the dimensionless Kerr 
parameter. Because of axial symmetry and stationarity, the Hamilton-Jacobi equation (1) 
for the Kerr metric has two cyclic coordinates, t and 0; therefore, Eq. (1) with the first 
integrals p t = —E = const and = L = const for the Kerr metric can be written as 

(1 - E 2 )((r 2 + a 2 ) 2 ~ Aa 2 sm 2 6)+ (4) 
AarLEGM/c 3 + (p 2 - 2GMr/c 2 ) x 
L 2 /(c sin 6) 2 A 2 S 2 +AS 2 = 2r (r 2 + a 2 )GM/c 2 . 

The existence of the complete integral of the Hamilton-Jacobi equation established by 
Carter (1968) is less obvious: 

S = -Etc 2 + 14+ j ^/0(&)d9 + j y^W)^, 



Q(9) = Q - cos 2 9(a 2 c 2 (1 - E 2 ) + L 2 /(sin 9) 2 ), (5) 
R(r) = (cE(r 2 + a 2 ) - La) 2 - (6) 
A(r 2 + (L — aEc) 2 + Q. 

Instead of the constant Q, we introduce a constant s < tt/2, which has the meaning of 
the minimum angle 9 (turning point) reached on the bounded trajectory in question [in the 
notation of Wilkins (1972), the constant s = 6L; see also Shakura (1987)]. According to (5), 
the constant Q can be expressed in terms of s: 

Q = cos 2 s(a 2 c 2 (l - E 2 ) + L 2 /(sin sf ). (7) 

In that case, the angle 9 on a bounded trajectory varies over the range s < 9 < tt — s. 

For stable trajectories on the r = const surfaces, the following conditions must be satis- 
fied: 

R(r) = 0, dR(r)/dr = (8) 
(Bardeen et al. 1972; Wilkins 1972). 



The newtonian analog of the Kerr solution 

The newtonian analog of the supercritical Kerr solution and its 
potentials 

In the Newtonian limit, E 2 ps 1 + 2H/c 2 and A «r 2 + a 2 . Substituting these expressions in 
Eq. (4) yields an approximate Hamilton-Jacobi equation: 

L 2 r 2 + a 2 

2H = (r 2 + a 2 ) sin 2 9 + r 2 + a 2 cos 2 9 S,r + (9) 

1 2 AarLGM/c 

]b,e + 



r 2 + a 2 cos 2 9 ' (r 2 + a 2 ) (r 2 + a 2 cos 2 9) 
2GMr A(GMrasm9) 2 



r 2 + a 2 cos 2 9 (r 2 + a 2 cos 2 9)c 2 (r 2 + a 2 ) 
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Note that Eq. (8) describes the trajectories of test particles in the field of a flat disk with 
radius a in special curvilinear coordinates r, 9, 0, which are related to Cartesian coordinates 

x, y, z by 

r = ^x 2 + y 2 + (z-ia) 2 + (10) 
+ y/x 2 + y 2 + (z + ia) 2 ), 

a sin 9 = — (\J x 2 + y 2 + (z + id) 2 — (11) 

\/ x 2 + y 2 + (z — ia) 2 ), = arctany/rr. 

When passing from Cartesian coordinates x, y, z to orthogonal curvilinear coordinates r, 9, <p 
using formulas (10) and (11), the nonzero metric tensor components are expressed in terms 
of the curvilinear coordinates as 

r 2 + a 2 cos 2 9 2 2 2 

^"■ = ' £e0 = r J + a 2 cos 2 fl, 12 

r 2 + cr 

= (r 2 + a 2 ) sin 2 0. 



Hamiltonian (9) contains the contravariant metric tensor components (12), and it can be 
written in Cartesian coordinates as 

= (p + A/c) 2 _ 
2 

The corresponding Hamilton equations are 

+ ^^ Pt + A t/ , (13) 

Eqs. (13) can be rewritten as the equations of motion in a gravitomagnetic field in the 
quasi-Newtonian approximation: 

dv 

— = V$ + curlA x v/c, (14) 

(JjL 

where the Newtonian potential $ and the gravitomagnetic vector potential A in a vacuum 
satisfy the Laplace equation. The equations of relative motion in a rotating coordinate 
system can be derived from Eqs. (14) if ft x R, R = (x,y,z), is substituted for A. In that 
case, curlA/c = — 2f2 and the second term on the right-hand side of (13) represents the 
Coriolis force. 

For a flat Kerr disk, the Newtonian potential $ is (Israel 1970; Zaripov et al. 1995) 

GM 

$ = — . + (15) 

2^x 2 + y 2 + (z + ia) 2 

GM 



2y^x 2 + y 2 + (z — ia) 2 



The vector A consists of the components —go a /goo, ot = 1, 2, 3 (Landau and Lifshitz 1980). 
In the Newtonian limit in a vacuum, the vector of the gravitomagnetic field \& can be 
introduced instead of the vector potential A: 

curlA = 2W. (16) 



The scalar \P in an axisymmetric case matches the imaginary part of the Ernst complex 
potential. For the special case of a Kerr disk, the scalar \1/ is (Zaripov et al. 1995): 



* = - (17) 

2i^/x 2 + y 2 + (z — ia) 2 

GM 



2i^/x 2 + y 2 + (z + ia) 2 ' 

Clearly, expressions (15) and (17) for $ and ^ satisfy the Laplace equation. 

As follows from definition (10), the r = const surface is neither a sphere nor an ellipsoid. 

The precession of circular keplerian Orbits in the field of a gravi- 
tating rotating mass 

For the Newtonian analog of a supercritical Kerr disk, a > GM/c 2 . However, since 
a < GM/c 2 for black holes, in this subsection, we discard terms of the order of a 2 in ex- 
pressions (15) and (17) for $ and \l/ lest the order of accuracy of the model be exceeded. In 
that case, $ « GM/r, * « -azGM/r 3 . 

Let us introduce a coordinate system with the origin at the gravitating mass with the 
z axis directed along the spin axis. The equations of motion for free particles in this coor- 
dinate system are 

du GMx ( , , z a\ GM , . 

dv GMy / . a\GM 

dw GMz ( . ^z \ GM . . 

Tt = ~— +(««(™-v = ») ? )— ■ (20) 

In our approximation, the complete integral of the Hamilton-Jacobi equation is [cf. for- 
mulas (5)-(6)] 

S = -Ht + L(j) + J ^Q{6)dd + J y f 2R{r)dr, 



Q(9) = L 2 cot 2 s- L 2 cot 2 9. (21) 
GM L 2 ^LaGM 
2r 2 sin 2 s r 3 c 



R (r) = H+ — - ^— -j- - 2—-. (22) 



In formulas (21) and (22), we chose the turning point for the angle 9 as a constant: according 
to (20), 9 varies over the range s < 9 < n — s, with < s < ir/2. 

For stable circular trajectories, the equalities R(r) = and dR(r)/dr = must be 
satisfied. Substituting expression (21) for a fixed radius r in these equations yields the 
corresponding energy H and angular momentum L: 

fl!s ™_ flv ^- (23) 

L w sin sVGMR - 3aGM^^. 

cr 



Consider the precession of circular orbits. From the Hamilton equations 9 = dH/ dpg and = 
= dH/dL, we have 

4- = 4\/cot 2 s - cot 2 9, (24) 
at r 2 

d<f> L 2aGM 



dt r 2 sin 9 r 3 c 



Integrating Eq. (24) yields 



/ 



Ldt 2aGMt . . 

+ —J— = (26) 



r 2 sin 2 9 r 3 c 



d9 2aGMt 
+ 



sin 2 9\/c 2 s - c 2 9 r 3 c 
cot A 2aGMt 



arcsm 



cot s J r 3 c 



When changing the variable in formula (26), we used Eq. (24). 

Consider the case where all particle at t — were in the sin sin 9 = cos 9 cot s plane or, 
introducing Cartesian coordinates, in the y = zcot s plane. According to equality (26), the 
particles will lie on the ycos (2aGMt/r 3 c) — a; sin (2aGMt/r 3 c) = zcot s surface at time t. 
The closer the particles to the rotating gravitating center, the faster their precession. For a 
given radius r, the particle motion in an accretion disk can be represented as the rotation of 
a circumference inclined to the z spin axis at angle s with angular velocity 2aGM/r 3 c. In a 
coordinate system corotating with the circumference, the particle rotates with a Keplerian 
velocity. Thus, if viscous friction between adjacent orbits is disregarded, a tilted accretion 
disk cannot be in a steady state. 



The energy and angular momentum of particles moving 
along geodesies orbits 

in the Kerr solution on the r = const surfaces 

In this section, for simplicity, we choose a system of units in which c — G — M — 1. To 

find stable orbits on the r = const surfaces, we must solve the algebraic equations (8). The 

corresponding solutions in the s = n/2 equatorial plane were found by Ruffini and Wheeler 

(1970) for j — 1 and by Bardeen et al. (1972) for arbitrary j. For s = and arbitrary j, 

the particle angular momentum L is zero. The corresponding expression for the energy was 

derived by Lightman et al. (1975). Inclined orbits at j = 1 were investigated by Wilkins 

(1972). In the case of charged rotating black holes the selected tilted bound geodesies were 

studied by Johnston and Ruffini in extreme case j — 1. In the case of charged rotating black 

holes the selected tilted bound geodesies were studied by Johnston and Ruffini in extreme case 

j = 1. Finally, Shakura (1987) derived formulas for the energy and angular momentum at 

arbitrary j and s by a complex, indirect method. Below, we show that derivation expressions 

for E and L from Eqs. (8) can be reduced to solving a quadratic equation and factorizing 

the numerator and denominator in the resulting fraction (to be subsequently canceled by a 

common factor). 

Let us introduce new unknowns instead of E and L: 

E E 2 — 1 

x c = T , y c = r2 ■ (27) 
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Eqs. (8) can then be written as 

y c a x + x c a 2 + (x 2 c - y c ){ai + a 3 ) + a = 0, 



y c b x + x c b 2 + (x c - y c )(bi + b 3 ) +b = 0, 
Peg 
dr 



bi — a \j,si * — 0,1, 2, 3. 



Here, 



ao = — (V 2 — 2r + A cos 2 s), 
ai = r 4 + 2j 2 + fr + j 2 A cos 2 s, 

a 2 = — 4jr, a 3 = — pA, p = r 2 + j 2 cos 2 s. 

Eliminating y from Eqs. (28) and (29) yields the quadratic equation 

xl{aib 3 - a 3 bi) + x c (a 2 b 3 - a 3 b 2 ) + 
+ a b 3 - a 3 b = 0, 

whose solution is given by 

E A + qB , ^~ 

x c = — = ——. , q = - j 2 cos 2 s r. 

L Dsms 

In formula (31), we use the notation 

A = j sin s(3r 4 - 4r 3 + fr 2 cos 2 s(j 2 r 2 - r 4 )), 
B = rpA, 
D = -a x b 3 + a 3 bi = (r 2 - j 2 cos 2 s) x 

/ 2 i -2\2 a -2 3 • 2 

(r + j ) — 4j r sm s. 

The numerator and denominator in (31) can be factorized as 

A + qB = r((r 2 + j 2 )q + 2jr sin s) x 
(j sin sq + r 2 + j 2 cos 2 s — 2r), 
D = r((r 2 + j 2 )g + 2jr sin s)((r 2 + j 2 )q — 2jr sin s). 

Therefore, cancelling a common factor in the numerator and denominator in (31), we 

E jq sin s + r 2 + j 2 cos 2 s — 2r 
c L (r 2 + j 2 )g — 2jr sin s) sin s 

Multiplying Eq. (28) by &i and Eq. (29) by a x and subtracting (29) from (28) yields 

2 1 x(a 1 b 2 - a 2 &i) + a\bo - a bi 

X c-Vc-J 2 - ^ 

_ A 1 + ggi 
D sin s ' 



A 1 = 2jrp{2r A - 3r 3 + j 2 r 2 + cos 2 sj 2 {r - j 2 )), 
B 1 = p 2 (r 3 - 3r 2 + j 2 + j 2 r). 

The expression A x + qB x can also be factorized as 

A l + qB 1 = p((r 2 + j 2 )q + 2jr sin s) x (38) 
(2rjqsms + r 3 — 3r 2 + j 2 cos 2 s(r + 1)). 

Hence, 



D sin s 



(39) 



_ A 1 + qB 1 
{ij' 2 + j 2 )q — 2jr sin s) sin s 
a/p(p — 3r + j 2 cos 2 s/r + 2jgr sin s) 

Using formula (36) for E/L, we obtain with (39) 

E = P-2r + jqsins (4Q) 

ypip — 3r + j 2 cos 2 s/r + 2jgr sin s) 

Formulas (39) and (40) are the sought-for expressions for the particle energy and angular 
momentum in inclined circular orbits. 

We determine the coordinate radius of the marginally stable orbit, r*, from the condition 
of E given by (40) being at a minimum. Representing r* as a Taylor expansion in terms of 
j, we obtain r* as a function of j, s: 

r»si6- s i n s + •• ~ 6 — 3.266j sin s — (41) 

j 2 (0.5 + 0.1111 cos 2s) + j 3 (-0.2532sins - 
0.0567 sin 3s) + /(-0.1 196 + 0.0206 cos 2s + 
0.0162 cos4s) + j 5 (— 0.11 sin s - 0.02 sin3s + 
0.0025 sin 5s)... 

At s = 0, the first expansion terms for r* are 

r* w 6 - 0.6111j 2 - 0.0828/ - 0.211j 6 - 
0.0066j 8 - 0.0023j 10 + . . . 

Substituting the Taylor expansion (41) in (39) and (40) yields Taylor expansions for the 
binding energy and angular momentum in the marginally stable orbit: 

1 - E* w 0.0572 + 0.0321j sin s + j 2 (0.0131 - 
0.0087 cos 2s) + j 3 (0.0143 sin s - 0.002 sin 3s) + 
j 4 (0.0065 - 0.006 cos 2s + 0.0003 cos 4s) + 
j 5 (0.0087sins - 0.002 sin 3s) . . . , 

w 3.4641 - 0.9428?' sins- j 2 (0.1123+ (42) 

sin s 

0.1443 sin 2 s) - j 3 (0.1178 sins + 0.0175 sin 3 s) + 

j 4 (-0.0138 - 0.0878 sin 2 s + 0.0138 sin 4 s) + 
j 5 (-0.0327 sin s + 0.0467 sin 3 s + 0.016 sin 5 s) . . . 
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We supplement expansions (41) and (42) with the values of the corresponding functions 
at j — 1 numerically constructed by least squares for an arbitrary fixed orbital inclination 
to the black-hole spin axis: 

1 - E* « 0.4222 - 0.8314 cos 2 s + 0.264 cos 4 s + 
0.8605 cos 6 s - 0.656 cos 8 s, 

r* « 1 + 6.2005 cos 4 s - 
7.1816 cos 8 s + 5.1365 cos 12 s, 

L* ~ sin s (1.1876 + 2.5913 cos 2 s — (43) 
0.9948 cos 4 s + 0.494 cos 6 s) . 

Note that for s = and j = 1, the radius of the marginally stable orbit is 1 + \/3 + 
^3 + 2^ « 5.2745 (in units of GM/c 2 ). 

Periastron and nodal precession of orbits on the r = 
const surfaces 

Choose the angle 6 as a parameter on a bounded trajectory. It then follows from (5)-(7) 
that 

d4 _ (2jEr + L(A/sm 2 6 - j 2 )) 

e = (l-^) (sin 2 ^ 2 (l-E 2 ) + -^), 
sin \ sin s J 

dt -2jrL + E(r 2 + j 2 cos 2 9) A + 2rE(r 2 + j 2 )) 

Expressions (39) and (40) for the energy and angular momentum in circular orbits must be 
substituted in these formulas for E and L. 

Integrating the right-hand parts of expressions (44) over 9 from s to it/ 2 yields an expres- 
sion for a quarter of the variation in azimuthal angle A0 and for a quarter of the period T in 
which the particle runs from a minimum latitudinal angle s to its maximum ir — s and back. 
The integration result can be expressed in terms of the elliptic functions that contain E, L, 
and s as parameters [instead of the constant Q, we inserted the constant s into formula (7), 
which has the meaning of orbital inclination to the black- hole spin axis). Using the notation 
for the complete elliptic integrals of the first, K(x), second, E(x), and third, n(n, x), kinds, 
we have 

T = AA (K(k) (jL - j 2 E + (Er 2 + Ej 2 - jL) x 
A0 = AA{m{- cos 2 s, k) + (2rEj - j 2 L)K(k)), 



in 



A={L 2 /sm 2 s + f{l-E 2 ))-^ 2 ; (45) 
k 2 = cos 2 sf(l-E 2 )A 2 . 

In contrast to Johnston and Ruffini (1974), who first wrote corresponding formulas in form of 
elliptic integrates for any j, in formulas (45), we inserted explicit expressions for all quantities 
in terms of the turning-point angle = s in the orbit. Expressions (39) and (40) must be 
substituted for E and L, respectively. 

For low-eccentricity orbitsQ(5r = esin£, e << r), it follows from the Hamilton-Jacobi 
equation (4) that 



d£ y/-Rf'/2 „„_d 2 R 

R = 7^ = ( 46 ) 



dd y/Q ' dr 

-2 ( 6(1 - E 2 )r 2 - 6r + (1 - E 2 )f sin s 2 + - 



L 



2 



sin 2 s 



From Eq. (46) for the periastron rotation frequency v r , we have 



u r = 2nu r = ^ = AA^^K(k). (47) 

Note that the frequencies A0/27rT and 1/T are commonly denoted by v§ and u$, respectively, 
so the sought-for nodal precession frequency is z/ no d = v$ — vq. The periastron precession 
frequency of an orbit is z/ per = u$ — v r (Merloni et al. 1999). 

To clearly present the result, let us derive asymptotic formulas from (45) and (47) in the 
form of Taylor expansions in powers of j by taking into account the dependence of the first 
integrals E, L, and Q on j and s given by formulas (39), (40), and (7). 

The constants E and L enter into the formulas for A0 and T via the ratios E/L and 
(1 — E 2 )/L 2 , Q/L 2 . Let us write out the series expansions of these ratios with the accuracy 
that will be required to calculate the nodal precession frequency (A<J)/2tt — 1)/T up to terms 
of the order of j 3 inclusive: 

E r-2 3r-4 



L r 3 / 2 sin s ~^ ^ r 3 

.Jr + 0.5r 2 + sin 2 s(8 
3 



r + 0.5r 2 + sin 2 s(8 + 7r - 1.5r 2 )) 



r 9/2 g j n s 

-E 2 r-4 ,8(r-2) 
L 2 r 3 sin 2 s ^r 9 / 2 sins 



Q 9 / 9(^ — 4) sin s , , 

§ = cot 2 ,(l +J 2 ^l (48) 

,4(r - 2) sins . 



r 



4 Syer and Clarke (1992) considered the possible existence of stationary disks in the equatorial plane 
in which the particle orbits were constant-eccentricity ellipses. Clearly, allowance for general-relativity 
effects will result in the intersection of orbital trajectories in such disks, because the orbital periastron 
precesses. However, in the Newtonian theory, unstable disturbance modes with a radially variable eccentricity 
disturbance exist in such disks (Lyubarskij et al. 1994). Arbitrary accretion-disk disturbances apparently 
produce spiral waves in the disk structure (Spruit 1987). 



I I 



Using (48), we obtain from expressions (45) and (47) for A0, T, and A£ 

27- 1 = rVJ-2^ 3(r - 4)smS + 
2-1.5r + (18-7.5r)sin 2 s) + ... 



,9/2 



(8 + 3r + (24 - 9r) sin 2 s) + . . . 
AA 2 r-6 , 12j, .„ . j 2 



+ — ^r(r - 2)sins^- x (49) 

3 15 \ 

-r 2 + 15r - 12 - sin 2 s(r - 4)( — r - 21) ) + . . . 

2 2 y 

In the limit s — > 0, when the spin axis becomes parallel to the tangential plane to the 
accretion-disk surface, the following formulas hold: 

^-l-^- + il(2-15r) + 

3j 5 /3 7 19 2 N 
1 - -r H r 2 I + 



r i5/2 ^2 4 32 
1^ + ^(8 + 30 + 



(^) ^ + ^( r2+10r - 8 > + < 50 > 
' -352 + 576r - 286r 2 + 9r 3 ) ] + . . . 



32r 7 



From formulas (49) for the Keplerian frequency u^, the nodal precession frequency u noA) and 
the periastron rotation frequency v r of a low-eccentricity orbit around the black-hole spin 
direction, we derive the compact formulas 

A0 l( (2 - 3 sin s)j 
U4 = 2ttz/ = _ = _ (^1 + _ + 

((12 + 9r) sin 2 s - 8 - 9r)j 2 



4 r 6 



+ . . . 



2j 2 1.5 sins 
u: nod = 2nv nod = (A0 - 2tt)/T w — - j — h 

,-2-3r + (6 + 1.5r) sin 2 s 

f Te + ••• 



M 2 = (27rz/ r ) 2 = r —^ + sin s(r + 2) + (51) 



2 



^(12 + 20r - 3 sin 2 s(l + r)(10 + r)) + ... 
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The first term in the formula for nodal precession was found by Lense and Thirring (1918) (see 
also Wilkins 1972). In a coordinate system rotating with the frequency v no d, the Keplerian 
orbit is stationary. 

In the limit s — > 0, when the black-hole spin axis becomes tangential to the disk surface, 
we can derive the following expansions for the nodal precession frequency, the periastron 
rotation frequency, and the Keplerian latitudinal frequency from formulas (45) and (47) 
using (50)Q 

2-nod = (A0 - 2n)/T « 2 1 (l - + 

, 4 8 + 54r + 27r 2 
3 16V6 + ' ' 



v — 6 
(2^) 2 = — + 



4-r 3?' 4 

-4" - -^(r 3 + 128r 2 + 172r + 32) + . . . 



2 ^ = ^-^( 3 - + 8 )+ (52) 
J '16 + 336r + 15r 2 ) + . . . 



g4 r 15/2 



To compare formulas (51) and (52) with observations, we should eliminate the radius from 
them and express the frequencies z/ no( j and v r as functions of the Keplerian frequency The 
functions z/ nod (z/^) and , v r {y^) may be said to be given by (51) and (52) in parametric form. 
An increase in z/ 2 and a decrease in the frequency difference between the two peaks z/ 2 — v\ 
during observations implies the transition of a radiating clump to an orbit closer to the black 
hole and the approach of its orbit to the marginally stable orbit. When observations are 
accumulated, formulas (51) and (52) make it possible to determine the orbital inclination 
to the spin axis of a slowly rotating black hole or neutron star. When analyzing the energy 
release in the boundary layer and in an extended disk for a rapidly rotating neutron star, we 
should took into account the appearance of an intrinsic quadrupole moment that exceeds the 
Kerr one (see Sibgatullin and Sunyaev 1998, 2000a, 2000b). We emphasize that formulas (51) 
and (52) are not related to the weak-field approximation and are valid up to the marginally 
stable orbit. 

If the disk almost lies in the equatorial plane of a black hole (s = vr/2), then formulas (51) 
give Taylor expansions of the formulas by Okazaki et al. (1987) and Kato (1990) (see also 
Merloni et al. 1999): 

1 

^4, 



7™3/2 _|_ j 



1 - v/1 - 4j/r 3 / 2 + 3j 2 /r 2 

Unod — 



.3/2 



+ j 



, 1 - 6/r + 8j/r 3 / 2 - 3j 2 /r 2 



W r- ( r 3/2 + j)2 • ( 53 ) 



5 Caution must be exercised when passing to the limit s — > 0; for example, when s — > 0, L/(yQ sin 2 I 
S(sm6), where 6(x) denotes the Dirac 6 function. 
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The nodal frequency of the marginally stable orbit 
for an arbitrary orbital inclination angle s 

The following relations (Sibgatullin and Sunyaev 1998, 2000a) hold for the marginally stable 
orbit in the equatorial plane in a Kerr field, in which the reciprocal radius of the marginally 
stable orbit acts as a parameter: 



. 4 Jx - V3 - 2x ^ /3 - 2x 

j = — & — ' E= V^' 

2 

(2y/xV3-2x + x), x = l/r. 



3\/3x 

Let us substitute these expressions in (53): 

_ 3x^ 2 
U(t> ~ 3- v ^v / 3-2x + 4x' 

uJnod = ^(1 + \ ~ V3 - 2x)), v r = 0. (53a) 

If the parameter x is eliminated from (53a), then we derive the following dependence by 
least squares: 

iv nod « 0.1872 - V6/36) + 
1.7246(w - v^/SG) 2 + 1.2064(ccv, - VQ/3Q) 3 . (53b) 

We emphasize that formula (53b) accurately describes the limiting dependence of the nodal 
precession frequency on Keplerian frequency at s = tt/2 in the marginally stable orbit over 
the entire range from 1/26 to 0.5, which corresponds to the range of Kerr parameter 
from —1 to 1. To derive dimensional dependences, we must make the following substitutions 
in all formulas: 23.1u; nod = u nod (M/lAM Q ), 23. lu;^ = i^(M/1.4M ), . . . , where frequencies 
v 4>i ^nod w iH be given hereafter in kHz. 

To determine the nodal precession frequency of the marginally stable orbit for an arbitrary 
orbital inclination s, we use formulas (45) in which for E, L, r, we substitute their expressions 
in the marginally stable orbit as functions of j and s, (41) and (42). Below, we give the highly 
accurate approximation dependences for the nodal precession frequency and the Keplerian 
frequency of the marginally stable orbit, in kHz, on Kerr parameter in the entire j range 
(— 1 < j < 1) that we derived by analyzing our numerical calculations for various inclinations 
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of this orbit to the black-hole spin axis: 



' = : Kod ~ ^^(0.2139j + 0.0448J 3 + 

0.0308j 5 + 0.007j 7 ), 
v+ « L ^f° (1.5716 + 0.1933j + 0.1539j 2 + 
0.0712j 3 + 0.0852j 4 + 0.0533j 5 ); 

n 1.4M . 9 

= 10° : < od ~ -^(0.2139j + 0.0431J 2 + 

0.0674j 3 + 0.054j 4 + 0.0004j 5 - 0.0414j 6 + 

0.038j 7 + 0.0432j 8 ), 

1.4M . „ 
v+ « ° (1.5716 + 0.3875j + 0.2206j 2 + 

0.0675j 3 + 0.1681?' 4 + 0.1289j 5 ); 

o * _ l-4M , noion , , 0.1394j% 
s = 30 : i/ wod « -^^(0.2139j + ^ _ - — ), 

* 1 4M 0/noion , 0.1923j% 
s = 45 : i/ wod « (0.2139j + - _ — ), 



* 


1.4M 


V nod ~ 


M 




1.4M 


^nod ~ 


M 




1.4M 


7 /* rvj 


M 


7 /* 


1.4M 




M 



60°: I ^«^(0.2130j + ^g) > 

80°: v* » (0.2139j + °- 2535 ^ ), (54) 

nod M \ J 1 _ o.976j' ; V ; 

Note that, in contrast to the preceding formulas, the last two formulas in (54) approximate 
the numerical data in the range — 1 < j < 0.99. 

The corresponding limiting dependence of the nodal precession frequency on Keplerian 
frequency for s = 7r/6, 7r/4, tt/3, 4ti/9 can be derived from numerical calculations by least 
squares (here, oJ^ = oj^ — \/6/36): 

s = 30° : w nod S3 0.3189cJ^ + 0.5205(cJ^) 2 + 

4.37(cJ ) 3 , 0.049 < ujj, < 0.171; 
s = 45° : ccw ~ O.2416cJ + 1.5953(cJ ) 2 - 

0.2723(cJ^) 3 , 0.044 < cj < 0.257; 
s = 60° : ~ 0.2031^ + 1.7744(cJ^) 2 + 

0.5149(^) 3 , 0.041 <uj^< 0.414; 
s = 80° : ccw ~ 0.1846cJ^ + 1.7929(cJ^) 2 + 

1.0862(cJ^) 3 , 0.039 < <= 0.485. (54a) 

In formulas (54a), we give the ranges of Keplerian frequencies that correspond to the range 
of Kerr parameter from —1 to +1. Consider an example. Let the black-hole mass be 
2.2M and the Keplerian rotation frequency be 1.2 kHz. Since the rotation frequency in the 
marginally stable Keplerian orbit of a nonrotating black hole with M = 2.2M is 23.1 x 
1.4/2.2 x v^/36 « 1 kHz, = 0.2 x 2.2/(1.4 x 23.1) « 0.0136. According to (54a), the 
nodal precession frequency of the marginally stable orbit is: z/ no d ~ 42 Hz for s = 7r/2; 
Vnod ~ 52.6 Hz for s = 7r/4; z/ no d ~ 65.3 Hz for s = 7r/6; z/ no d ~ 93.6 Hz for s = vr/18; and 
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^nod ~ 123 Hz for s — > 0. We see that the nodal precession frequency significantly depends 
on the orbital inclination and changes by a factor of 3 as the inclination changes! If the nodal 
precession frequency is identified with the horizontal-branch oscillation frequency ^hbo> then 
it becomes possible to determine the inclination of the marginally stable orbit from the known 
black-hole mass, the Keplerian frequency, and the nodal precession frequency At j << 1, 
the nodal precession frequency of the marginally stable orbit is related to the Keplerian 
frequency by 

2 

uJnod ~ r , » . — v<j>. (54b) 
2 + 9 sin s 

At a fixed inclination s, the precession frequency v at the marginally stable orbit mono- 
tonically increases with j and reaches a maximum at j — 1. For the mass M = 1.4M and 
s = 0, it is 296 Hz. 

The larger the inclination s at fixed j, the larger the precession frequency f no d- As 
an illustration, let us write out the approximation formula for z/ no d at j = 0.998 [Thorne 
(1974)fl 

17.1094 sin 6 s , . , 

v nod « 0.29 + 0.7783 sin s + ^in kHz. 55 

1 + 1.35 sin 4 s V ; 



Conclusions 

According to the ideas explicitly formulated by Bardeen and Petterson (1975) and Bardeen 
(1977), a geometrically tilted accretion disk is most commonly modeled as a set of rings 
with the center at the coordinate origin that smoothly turn with changing radius under the 
effect of viscous torques and gravitomagnetic forces. This approach underlies the existing 
theories of tilted accretion disks [Papaloizou and Pringle (1983) and Markovic and Lamb 
(1998, 2000) for small disk deviations from the equatorial plane; Pringle (1992, 1997) for 
finite disk deviations from the equatorial plane), which are based on a vector equation for the 
conservation of angular momentum. Some ideas developed by Shakura and Sunyaev (1973, 
1976) for flat disks can be used to study the physics of tilted accretion disks. Van Kerkwijk 
et al. (1998) explained the puzzling spindown and spinup of some X-ray pulsars by the fact 
that the accretion-disk tilt in the inner regions could become larger than 90 degrees! 

The effects considered above take place for the inner nonstationary part of the disk, where 
viscous torques may be disregarded. Matter is not accumulated in the marginally stable orbit 
because of its instability. In contrast to the statement by Stella et al. (1999) that the orbital 
inclination affects weakly the nodal precession frequency, our detailed analysis shows that 
this effect is significant. There are probably preferential inclinations of the marginally stable 
orbits, and there is no need to introduce even harmonics of the nodal frequency to reconcile 
observations with theory. 

6 Studying the evolution of a black hole under the effect of disk accretion in the equatorial plane, Thorne 
adduced arguments that the Kerr parameter cannot be larger than 0.998 in this case when the radiation from 
inner disk regions is taken into account. Since the black hole mostly captures photons with a negative angular 
momentum (opposite to the black- hole rotation), the black- hole spinup by disk accretion is counteracted by 
its spindown via the predominant capture of the photons propagating against the rotation. 
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